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60 STONE. ON THE COMPUTATION OF THE POSITION IN ORBIT 

ON THE COMPUTATION OF THE POSITION IN ORBIT OF A PARTICLE IN 

THE TAIL OF A COMET. 

By Prof. Ormond Stone, University of Virginia. 

We have the well known equations 

''■~*(H]' (i) 

k 

f, cos /?, = — — e sin v, (2) 

/,sin£,=^; (3) 

where, neglecting the mass of the comet, 

k is the constant of the solar system, 

e, a, p, the eccentricityysemi major axis, and semi parameter of the orbit of 
the comet, 

r h v h the radius vector and true anomaly for any given time t,, 

/,, the corresponding velocity in orbit, and 

,9,, the angle which the radius vector makes with the tangent. 

/ \ j t \ ■ «. a P l x + e cos v < 1 \ 

(2) and (3) give tan B,— i -. — : = : . (a) 

K ' w/ " r r, e sin v, e sin v, v ' 

For a free particle moving with the comet up to the time t,, but after that 
moving in a hyperbolic orbit subject to a repulsive force Px acting in addition 
to gravity, but in opposition thereto, we have 

/■=*»(,_*) (!+}], (4) 

fco a p=* y^- x) Esmr, (5) 

f sin /9 = r > ( 6 ) 

in which /; is the hyperbolic semi-major axis. 

(1) and (4) give ? - l - = (1 - x) [ 2 - + | J , 
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since at the time t,, f =f, and R = r, ; 



I — x r, 



.-.*=- -.— iJ— , (J) 

/z I -%wc v y 

I X 

= — — r„ when the orbit of the comet is a parabola. 
Putting £ = sec <p, whence P = <$tan 2 ^, (i) and (6) give 

/,sin/3, = &^ [^ — ij . sin/?, = £,/(i — x)J^-tan^; 
whence, substituting the value of b given in (b) and solving, we obtain 

2i/x 

= — - — sin/?,, when the orbit of the comet is a parabola. 

The sign of E is assumed the same as that of tan <p, hence negative when x > i. 

(5) and (6) give tan fi = ^ . R -^y , ( 7 ) 

where, if F be the hyperbolic anomaly, 

R sin V = b tan <p tan F. 

Substituting in (7) the values of E, P, and A 1 sin V, and simplifying, we have 

tan F = sin <p cot /9 ; (</) 

or, for the time /„ tan F, = sin ^ cot /9,. (<^,) 

After computing the formulae (a), (b), (c), and (</,), r, and F, are found from 
the well known equations 

r, sin V, = b tan </> tan T* - ,, 
r, cos V, = b(E — sec F,). 

The agreement of this value of r,. with that derived directly from the orbit of 
the comet will form a convenient check. 

The mean motion and time of hyperbolic perihelion passage t is found 
from the equations 

kky'ii —x) 

N, = n{t, — t ) = XE tan F, — log tan (45 ° + %F,), 
where X is the modulus of common logarithms. 
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The position in orbit for any given time t is found from the formulae 

/ E tan F — log tan (45 ° + \F) — n (t — /„) = N, 

Rsin V = b tan <p tan F, (/) 

R cos V =b {E — see F). 

F is found from the first of these equations by successive approximation. 
If £sec F be negative or be positive and large, i/(EsecF — 1) will be small 
and we may obtain tan F x from the equation 

tan F - N + l ogtan(45°+^ ) ) 
tan t x — j £ 

where F is an assumed value of F. Subtracting this from 

tan F =K+A°z t*nj 4J1± ¥) 

AE 
we have, including only terms of the first order. 



F=F l + 



E sec F — 1 ' 

If F be not too large we may express the first of equations (e) in the form 
E . \ (e< - e~*) - x = », {t - t a ) = iV„, 
where Ax = log tan (45 + %F} 

and n = -?-^ '- . 

„. e* — e~* x , x 3 , jr* 

Since ="H r- • • •, 

2 1 I .2.3 I .2.3.4.5 

if we put 1/ = 1 -J 1 - (-..., 

we have after reduction 

r s . 6(^-1) J_ 67V 1 

- r + £ ■ fi x £ ~t?-°' 

a cubic equation, which may be solved by means of the following formulae: — 
u = 2 (E — i)/E, «, = E/3N0 ; 
a = «;,«,*, ,9 = 2«o* ; 

tan = «//i, tan f = tan* -|#, .r = ,?//./ tan f . 

/i is obtained from the accompanying table with the argument x 2 . In the first 
approximation, unless otherwise known, ;i may be assumed equal to unity. 



SOLUTIONS OF EXERCISES. 



63 



X 2 


log 11 


Diff. 


x 2 


log fl 


Diff. 


X* 


log /I 


Diff. 


X 2 


log /l 


Diff. 


0.0 

O.I 

0.2 

0.3 
0.4 


0.00000 
0.00109 
0.00217 
0.00326 

0.00434 


109 
108 
109 
108 


1.0 
1.1 
1.2 

'■3 
1.4 


0.01084 
0.01 193 
O.QI30I 
0.01409 
0.015 17 


109 
108 
108 
108 


2.0 
2.1 

2.2 
2-3 
2-4 


0.02166 
02274 
0.02382 
0.02490 
0.02598 


108 
108 
108 
108 
108 
108 
108 


30 

3-1 

3-2 

3-3 
3-4 


0.03245 
0.03353 
0.03460 
0.03568 
0.03676 


I08 
107 
108 
108 
107 
108 


0.5 


0.00543 


109 
108 
108 


••5 


0.01626 


108 
108 
108 
108 
108 


2-5 


02706 


3-5 


0.03783 


0.6 
0.7 


0.0065 1 
0.00759 


1.6 
i-7 


0.01734 

0.01842 


2.6 

2-7 


0.02814 
0.02922 


3-6 

3-7 


0.03891 
0.03998 


107 
108 


0.8 
0.9 


0.00868 
0.00976 


108 
108 


1.8 
1.9 


0.01950 
0.02058 


2.8 

2.9 


0.03029 

003137 


108 
108 


3-8 
3-9 


0.04106 
0.04213 


107 
108 


1.0 


0.01084 


2.0 


0.02166 


30 


0.03245 


4.0 


0.04321 
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The dome of the rotunda of the University of Virginia is spherical. The 
length of its meridian section is %s f .2 and the girth of its base is 2i4 / . From 
these data it is required to compute the radius and the surface of the dome. 

[William M. Thornton.] 

SOLUTION. 

We have for the chord of the meridian arc 

2 -i4 = 68.i. 

If therefore x be the half central angle of the arc, 
sin x 681 

_ _ _ 0.7QQ. 

x 852 '** 

Passing to logarithms and expressing x in minutes by putting 

x = x . 



10800' 

we get y = log sin- x — log x + 3.6336 = o. 

By trial we find, if 

x = 64°, y = + 0.0030 ; 



* = 65°, y 



0.0002 ; 



